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We present a detailed simulation of the performance of water Cerenkov detectors suitable for use in the Pierre 
Auger Observatory. Using GEANT4, a flexible object-oriented simulation program, including all known physics 
processes, has been developed. The program also allows interactive visualization, and can easily be modified for 
any experimental setup. 



1. Introduction 

Water Cerenkov detectors have proved to be 
superb devices for the study of cosmic air show- 
ers [Q, and will constitute major components of 
future experiments such as the Pierre Auger Ob- 
servatory H . Monte Carlo simulation of the re- 
sponse of such detectors to signal and background 
processes is crucial for proper interpretation of 
the data and to aid in development of reconstruc- 
tion and analysis algorithms. To address these 
needs, experiments have prepared dedicated de- 
tector simulations (^]. There may be advantages, 
however, in exploiting some of the efforts which 
the High Energy Physics (HEP) community has 
directed at this problem. The GEANT3 pack- 
age for example, was developed in order to 
provide HEP experiments with generic tools for 
simulating the passage of particles through mat- 
ter, but it also found substantial usage in the 
medical and biological sciences and in astronau- 
tics, and in fact has been studied for use in sim- 
ulating water Cerenkov detectors for the Auger 
experiment ||. In 1996 CERN initiated the 
GEANT4 project with the goal of reproducing 
all the functionality of GEANT3 using an Object 
Oriented approach, as well as addressing some of 
the shortfalls of the older program. We have be- 
gun investigating the suitability of GEANT4 for 
the problem of simulating water Cerenkov detec- 



*doqui@hepmail. physics. neu.edu 
tmccauley@hepmail.physics.neu.edu 
t torn, paul@hepmail.physics.neu.edu 
8 Stephen. reucroft@cern.ch 
^john.swain@cern.ch 
II lucas.taylor@cern.ch 



tors, and report here on the status of the work. 

2. GEANT4 

Like its predecessor, GEANT4 provides a bat- 
tery of tools to describe the geometry and ma- 
terial properties of an experimental setup, han- 
dle particle transport through materials and mag- 
netic fields, and simulate particle decay and in- 
teractions with detector elements. All relevant 
physics processes have been included. For ex- 
ample, simulation of an electron traversing a 
material can include effects of ionisation, delta 
ray production, multiple coulomb scattering, 
bremsstrahlung and Cerenkov radiation. Optical 
photons produced by processes such as Cerenkov 
radiation may then be subjected to Rayleigh scat- 
tering, absorption, and optical boundary interac- 
tions. We note that the calculation of reflection 
and transmission coefficients at material bound- 
aries takes into account the polarization state of 
the photon, which is important in accounting for 
a photon's fate as it traverses multiple optical 
boundaries en route to a photocathode. 

In addition to reproducing the functionality 
of GEANT3, GEANT4 aims to improve the proce- 
dures used for geometry definition, introduction 
of special physics processes, visualization, and op- 
tical processes. In addition, special effort has 
been made to ensure tracking precision over ar- 
bitrary scales, so the package may lend itself to 
simulation of extensive air showers |p"l|| as well 
as detailed detector response to ground particles. 
GEANT4 is written in C++ using an object ori- 
ented approach, and as such alleviates some of 



2 




the shortcomings inherent in the procedural ap- 
proach (FORTRAN77) of geant3. 

3. Detector Geometry 

GEANT4 strives to provide rather advanced 
tools for describing detector geometry, such as 
methods for interpreting files produced by Com- 
puter Aided Design (CAD) systems. However, 
for our application a first approximation to the 
geometry may be described in only a few lines of 
code. Figure [l] shows an example of an Auger- 
style Cerenkov detector which consists of a cylin- 
drical tank filled with water. A reflective liner 
surrounds the water volume. Three hemispherical 
domes at the top of the tank house photomulti- 
plicr tubes, which contain sensitive photocathodc 
volumes that register a hit when a photon strikes. 
Photons bounce about inside the tank until they 
are absorbed in the water, the tank walls, or un- 
til they enter a photomultiplier. The frequency- 
dependent optical properties of all these elements 
can be tuned to available data. In the figure, a 
1 GeV muon is incident at the top of the tank, 
and some of the radiated Cerenkov photons are 
shown. 



Figure 2. Distribution of reflected photons as 
a function of polar angles 9 and <f> for a beam 
incident at ^incident = 30°. This simulation uses 
the unified model of GEANT4 with specular and 
backscatter spikes, visible at <j> = 90° and 270° 
respectively, a specular lobe centered at = 30°, 
as well as uniform (Lambertian) reflection. 



4. Optical Modeling 

The detector geometry of Figure [l] includes a 
diffuse reflector to contain the Cerenkov light in- 
side the volume of water. The use of Tyvek [] as 
such a reflector has been explored by the Auger 
experiment, and it has been noted H that the 
treatment of photon interactions with a rough di- 
electric surface as available in geant3 does not 
appear to be sufficient to describe the experimen- 
tally observed reflection of light from this ma- 
terial. GEANT4 provides a more flexible optical 
model inspired by the work of Nayer et al. |||. 
There it was observed that the principal features 
of both physical and geometrical optical models 
of surface reflection could be accommodated in a 
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so-called unified model which is applicable over a 
wide range of surface roughness and wavelengths. 

The unified optical model in GEANT4 allows 
adjustment of parameters to control the relative 
contributions from: specular reflections about 
both the average surface normal and the nor- 
mal of a microfacet at the surface; the diffuse or 
Lambertian reflection; a possible backscatter con- 
stant; overall surface reflectivity. Figure |] shows 
the polar angular distribution of light bouncing 
off such a unified surface, with contributions from 
the various sources identified. Tuning of this 
model to experimental data j?]] is underway. 

5. Benchmarking 

On a 450 MHz Pentium, our program re- 
quires roughly 9 seconds to simulate the tank re- 
sponse to a 1 GeV vertical muon|. We anticipate 
substantial improvement in this figure, possibly 
through application of top-C j[o| or parameteri- 
zation tools built into the GEANT4 framework. 

6. Summary and Prospects 

A prototype simulated water Cerenkov detec- 
tor has been developed using geant4. We are 
optimistic about the utility of this package for a 
number of reasons. 

• The general, extensive arsenal of tools pro- 
vided will allow the detail inherent in the 
simulation to evolve as needs are identified. 
On the other hand, fast parameterizations 
can be employed in a straightforward man- 
ner wherever details might be found to be 
unnecessary. 

• GEANT4 design is consistent with modern 
software trends. Though the product is rel- 
atively new and potholes in the road should 
be expected, the project enjoys a large de- 
velopment community, which bodes well for 
its longevity and evolution. 

• Several shortcomings of GEANT3 have been 
addressed in the new version. For example 
the surface reflection model now available 



is more amenable to our problem than its 
predecessor. 
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